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Domain-walls in one-dimensional Ising ferromagnets can undergo Bloch oscillations when sub¬ 
jected to a skew magnetic field. Such oscillations imply finite temperature non-dispersive low- 
frequency peaks in the dynamical structure factor which can be probed in neutron scattering. We 
study in detail the spectral weight of these peaks. Using an analytical approach based on an approx¬ 
imate treatment of a gas of spin-cluster excitations we give an explicit expression for the momentum- 
and temperature-dependence of the spectral weights. Generally the spectral weights increase with 
temperature T and approaches the same order of magnitude as the spin-wave spectral weights at 
high temperatures. We compare the analytical expression to numerical exact diagonalizations and 
find that it can, without any adjustable parameters, account for the T and momentum-transfer 
dependence of the numerically obtained spectral weights in the parameter regime where the ratio of 
magnetic fields h x /h z <g 1 and the temperature is h x < T <~ J z / 2. We also carry out numerical 
calculations pertinent to the material CoNb 2 C> 6 , and find qualitatively similar results. 

PACS numbers: 


I. INTRODUCTION 

Bloch Oscillations (BO) 1 2 is a pure quantum mechan¬ 
ical phenomenon where a particle in a periodic poten¬ 
tial oscillates when acted upon by a constant external 
force. Although BO are quite delicate, they have been 
observed in a range of systems in the last two decades 3 6 . 
In addition BO have also been predicted- to exist in one 
dimensional magnetic systems, so called magnetic Bloch 
oscillations (MBO). However, such MBO have so far not 
been observed experimentally. 

One way to search for MBO is to search for its energy 
spectrum which is the Wannier-Zeeman (WZ) ladder^ 
of equally spaced energy levels. Experimentally one can 
hope to probe the WZ ladder directly using neutron scat¬ 
tering. At low temperatures neutrons will induce transi¬ 
tions from the ground state to the WZ levels. However, 
the transition matrix elements of these processes are very 
small, thus the low temperature inelastic neutron scat¬ 
tering signal will be very weak 1 ^ 11 . In contrast, at tem¬ 
peratures comparable to the dominant magnetic scale the 
WZ levels will be thermally populated and neutrons can 
induce scattering between them. These transitions show 
up as peaks at low frequencies in the neutron scatter¬ 
ing dynamical structure factor, and are the topic of this 
paper. 

This paper focuses especially on the spectral weights of 
the finite temperature low-frequency peaks and how they 
depend on momentum transfer and temperature. These 
quantitative details are important in order to distinguish 
the peaks associated with BO from generic low-frequency 
peaks in other systems with Zeeman spin-split spectra. 

In order to calculate the spectral weights we use both 
an analytical and a numerical approach. The analyt¬ 
ical approach is based on a gas of spin-cluster excita¬ 
tions with collisions treated in an approximate way. This 
goes beyond the one 7 ^ and two 10 domain-wall approxima¬ 
tions where collisions are neglected completely. The ne¬ 


glect of collisions can be expected to work well at low 
temperatures, but is more dubious at high temperatures 
where one expects frequent collisions between domain- 
walls. Also the one- and two domain-wall approximations 
restrict the state space severely to states that constitute 
the WZ levels. Since the existence of the low-frequency 
peaks is caused by the thermal population of the WZ 
levels it is important to verify that thermal population 
of other excluded states will not destroy them. For these 
reasons we validate our analytical result by carrying out 
numerical exact diagonalizations keeping all the states. 

While the analytical result presented in this paper is 
strictly valid for the Ising model in a skew magnetic field 
we expect it to hold qualitatively for a broader class of 
Hamiltonians. We substantiate that at the end of the 
paper by performing numerical exact diagonalizations of 
a Hamiltonian that also includes extra terms relevant for 
describing the material CoM^Og. 

The outline of the paper is the following. In chapter 
2 we present the Hamiltonian and explain qualitatively 
why it should describe MBO. Chapters 3 and 4 explain 
the analytical calculation of the dynamical structure fac¬ 
tor. The numerical exact diagonalization result is pre¬ 
sented in chapter 5, while chapter 6 is devoted to a de¬ 
tailed comparison of the analytical and numerical results. 
In chapter 7 we present numerical results for the dynam¬ 
ical structure factor for parameter values relevant to the 
material CoNb 20 g, and chapter 8 concludes the paper. 


II. HAMILTONIAN 

The prediction that domain-walls in one-dimensional 
ferromagnets undergo BO in a magnetic field is rather 
general and holds for a range of Hamiltonian^. In or¬ 
der to be definite we will focus on the one-dimensional 
spin-1/2 ferromagnetic Ising Hamiltonian with nearest- 
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neighbor couplings in a uniform skew magnetic field 

H = ~ E 1 + /fcSf + ■ (1) 

i 

We will consider the case where the Ising coupling J z 
is the dominant one, and J z h z > h x > 0. When¬ 
ever numerical energies, frequencies or magnetic fields 
are quoted, they are in units of J z . The lattice spacing 
and h are set to unity. 

The ground state of the Hamiltonian (|T|) resembles 
closely the ferromagnetic all spins up state as quan¬ 
tum fluctuations are small. Excitations can be classified 
according to how many domain-walls they contain. A 
domain-wall is an antiferromagnetic arrangement of two 
neighboring spins and costs an energy J z j 2. It is these 
domain-walls that undergo BO in a skew magnetic field. 
To see why, consider first the transverse field h x . It cou¬ 
ples to the spin operator S x that flips a spin. When 
J z is large the main effect of the transverse field is to 
flip a spin adjacent to a domain-wall which causes it to 
move, and will broaden the domain-wall excitations into 
a band. Thus the role of h x is to mimic a periodic poten¬ 
tial. The longitudinal magnetic field, h z > 0, plays the 
role of the external force as it will pull the domain-wall 
in the direction of the down-spins. The resulting BO of 
a single domain-wall in this system was recently studied 
in real-time numerical simulation^!. 

However, in contrast to other systems with BO, a sin¬ 
gle domain-wall cannot exist in isolation in the magnetic 
system. The competition between the longitudinal mag¬ 
netic field and the Ising coupling means that the system 
cannot have arbitrarily long spin-down domains. Instead 
new domain-walls will be created, and the excitations 
will be pairs of domain-walls, i.e. domain-wall/(anti)- 
domain-wall bound states. Nevertheless, the magnetic 
field still causes the domain-walls to oscillate, thus MBO 
in domain-wall ferromagnets resemble breather oscilla¬ 
tions. 

III. SPIN-CLUSTER EXCITATIONS 

At h x = 0 the excitations are domains of l consecutive 
down-spins surrounded by up-spins. For h x > 0 this is 
slightly altered, and the excitations will be superpositions 
of such domains with different numbers of down-spin^®. 
These excitations have been termed multiple-magnon 
bound stated^, but we will call them spin-cluster excita¬ 
tions for simplicity. A spin-cluster state can be described 
mathematically by labeling the state with l consecutive 
down-spins starting at site j as | j, l). Defining Fourier- 
transformed states as | p,l) = -l=e _ * p */ 2 JT \j, l), 
and writing energy eigenstates as superpositions of these, 

| p, n) = i Vv>ra(0I.P> 0> the energy eigenvalue equation 

in the sector with two domain-walls takes the form 

h x cos (p/2) [4>p, n (l + 1) + ip p ,n{l ~ 1)] = 

(Jz + lh z — E p , n ) ip Pj n(l) , l = 1,2,... (2) 


where we have defined ^p#o,n(0) = 0, and neglected all 
terms that change the number of domain-walls. This is 
expected to be a good approximation for large J z . 

Eq. |2]) resembles the recursion relation for Bessel func¬ 
tions and has the solution 

Ep } n = J z -f- h z Vp n (3) 

and 

= E n Ji—i/ pn (Zp) ( 4 ) 

where C n is a normalization constant, n € {1, 2,...} la¬ 
bels the energy levels and z p = xcos(p/2). We have 
introduced the dimensionless ratio x = 2h x /h z . To get 
a normalized wave function in the infinite chain limit we 
have set the coefficient in front of the increasing Neu¬ 
mann function solution A /_ u n ( z p ) to zero. The num¬ 
bers v n n are obtained by the other boundary condition, 
Eq. §' for l = P! 

J-*A*r) = »■ (5) 

For z p >~ 1, v pn will depend on p for the lowest values 
of n, thus the lowest lying spin-cluster excitations will 
be dispersive 10 . For z p <~ 1, v PjTl « n. Thus for x < 1, 
E Pt n « J z +h z n is independent of p and the energy levels 
are equally spaced with a separation h z . This is the WZL 
and we expect finite temperature neutron scattering to 
reveal transitions within the ladder with characteristic 
frequencies w = kh Zl where k G Z. The wave function 
describing the spin-cluster excitations in this WZL is a 
superposition of domains of different lengths /, each with 
weight ip p , n (l) = C n Ji_ n (z p ) for / > 0. 

IV. DYNAMICAL STRUCTURE FACTOR 

We will now use the spin-cluster excitations to calcu¬ 
late the spectral weight of the low frequency peaks ana¬ 
lytically. The transverse dynamical structure factor is 

1 r°° dt 

S xx ( q ,u) = -J2 J ^e^(i\e-^ H S x q (t)S x (0)\i) 

(6) 

where Z is the partition function and /3 is the inverse 
temperature. The spin operators in momentum space 
are S x = e~ lqr S x and have time dependence S x (f) = 
e lHt S x e~ lHt . The sum is taken over all energy eigen¬ 
states of the system. The time integral can be performed 
by inserting an extra sum over energy eigenstates 

S**(q,w) = ±Y,5(w-( E i - E i ))e-f ,E ' |0K|i)| 2 . (7) 

ij 

We will assume that the energy eigenstates can be writ¬ 
ten as states of almost independent spin-clusters. Thus 
we specify an energy eigenstate with d spin-clusters 
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as \ni,pi,n, 2 ,P 2 , ■ ■ ■ ,nd,Pd) where n r labels the inter¬ 
nal quantum number of spin-cluster r, and p r its mo¬ 
mentum. we assume that this state has energy E = 
J z d + h z 1 n r■ Thus we neglect any interaction en¬ 
ergies between different spin-clusters. Such many-spin- 
cluster states are orthogonal and complete for x = 0. 
For finite x we expect orthogonality and completeness to 
hold approximately when the density of spin-clusters is 
not too high. 

The internal quantum number n r is an integer, so the 
dynamical structure factor will have peaks at low energies 
when Ej — E t = kh z for k £ Z. To get the spectral 
weight S£ x (q) of the peak at w = kh z we integrate over 
frequencies in the vicinity of kh z 

St x {q) = ^E e '^l01-5 9 1i)| 2 k- Bi =^ (8) 

ij 

where the sum over states is restricted to states \ j) which 
energies differ from |i) by kh z . 

When the operator S x acts on the state |i) it gives 
a sum over states | i') which each differ from state \i) 
by having a single spin overturned. Many of these 
states do not contribute in Eq. as they have an 
almost zero overlap with states \j) having energy Ej 
that differ from Ei by kh z . This is in particular true 
when S x creates(destroys) two domain-walls. This in- 
creases(decreases) the energy by J z and will therefore 
not contribute to low energy peaks at w = kh z . The only 
states |*') that contribute are those which is the result 
of flipping a spin so as to displace an already existing 
domain-wall. When the spin-clusters are not too densely 
packed we can treat the spin-clusters independently. In 
that case the matrix element reduces to the matrix ele¬ 
ment between two single spin-cluster states as the other 
spin-clusters are unchanged. However, at the tempera¬ 
tures considered here, one also needs to go beyond the 
independent approximation and take into account that 
the movement of a domain-wall can cause it to collide 
with the neighboring spin-cluster. The result of such a 
collision is the disappearance of two domain-walls and 
a merger of two domains. The energy of such a final 
state differs by the initial state by roughly J z , and should 
therefore not be counted as a contribution to the low fre¬ 
quency peaks. We thus need to explicitly exclude final 
states | i') where a domain-wall is moved such as to merge 
two domains. We achieve that by including projections 
P(T which meaning is that the domain r in state |i) must 
be such that its left domain-wall can be moved one lattice 
spacing to the left without touching another domain-wall 
for there to be a contribution. We can therefore write the 
matrix element that contribute to peaks at kh z as 

di 

E It? Wl 2 = E I A n rPr Pr L +At Pr P^ R +Bn rPr +Bl pr \ 

3 r=l 

( 9 ) 


where di is the number of domains in state |i) and 

1 OO 

AnrPr. = 2 EW + t(' + m Pr ,n r (l)e i(ql - p ^ )/2 , 

1=1 

(10) 

OO 

Bn rPr = 2 EC+ 9 ,rv+^ ^ l)Vw r (0e i(9 ' +Pr)/2 , 

(11) 

are matrix elements between two single-cluster states, 
one with quantum numbers ( n r + k,p r + q) and its 
left(right) domain-wall displaced one lattice spacing to 
the left compared to the other state which have quantum 
numbers ( n r ,p r ). The right movement is described by 
their complex conjugate. 

By multiplying out the right hand side of Eq. © we 
get 

= + p rp 

i r =1 ^ 

x [l^n^pj 2 + ( A nrPr + A* nrPT )(B nr . Pr + B* rPr )] 

[An rPr + A^ rPr ] + ( B nrPr + B^ rPr ) 2 ^j 

(12) 

where we have set P^~ L P^ = P^~ L and P^P^ = P 

The first term, with P^, is proportional to the prob¬ 
ability to find a spin-cluster with quantum numbers 
(ni,pi) whose left domain-wall can be expanded to the 
left without touching the neighboring domain (periodic 
boundary conditions is assumed). This probability does 
not depend on the spin-cluster index r, thus the sum 
over r can be taken into account by replacing the prob¬ 
ability with the number of occurrences. We approxi¬ 
mate this quantity with the number of occurrences of 
n consecutive down spins surrounded by one up spin 
to the right and two up spins to the left, = 

J2i(Ui-iUiD i+1 ■ ■ ■ D i+n U i+n+x ), where U z = 1/2 + Sf 
(Di = 1/2 — Sf) is the projection operator onto up(down) 
spins at site i. The brackets denote both quantal and 
thermodynamical average. The two up-spins to the left 
allow the displacement of the left domain-wall one lattice 
spacing to the left without touching a neighboring do¬ 
main. For x = 0 this approximation is exact. For small 
x it is valid approximately as the wavefunction i/„ iP (/) is 
dominated by the term with l = n. Similarly the terms 
with P r “/ is approximated by by N^ n ^ 

and the terms without projections by Nj- n q. This results 
in 

S% x (q) = EE { [ 2 I ^l 2 + p + K P )(Bn P + B* np )] N nnt 

n p ^ 

+ [An p + A^ p ] IVpi-ntt + (B np + B^ p ) 2 N^ n ^j 


(13) 
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where we have used inversion symmetry = N^ n -f. 

At high temperatures, T > h x , we can approximate 
the N^ n - |-s by setting h x = 0 and use transfer matrice^. 
In the thermodynamic limit the number of n down-spins 
surrounded by m up-spins to the left and l up-spins to 
the right is 

A T ... t4 ... 4t ... T = NP t K^~ 1 WK!!r 1 WK l + 1 (14) 

where P-j- = 1/2 + sinh(/3/i z /2)[2(A — cosh(/3/i z /2))] _1 is 
the probability of finding a certain spin to be up. W = 
g-pjz/z/\ j s the domain wall factor and K± = e ±/3ft, */ 2 /A 
is associated with finding two adjacent spins that points 
up(+) (down(-)). The constant 

A = cosh (f3h z /2) + [sinh 2 (/3ft. z /2) + e _/3J *]^ (15) 

is the largest eigenvalue of the transfer matrix. Similarly 
one can find the number of domain-walls that separates 
exactly n up-spins and to down-spins. It is in the ther¬ 
modynamic limit 

A t ... n .. 4 = N^K^WK™- 1 (16) 


where 7 = e -0-W 2 [2(A — cosh(/3/i z /2))] _1 . 

The expressions for A np and B np , Eqs. (10 1 and 
involve an infinite sum over l. In order to evaluate tins 
we add and subtract terms with negative values of Z, see 
Appendix. Substituting the resulting expressions into 
Eq. (13) gives 


S% x {q) = \N m J^ k {w q ) (Jr-fcK) - 2J_ 1 _ fc (u>,)) 

+ -k( w q) 

- \ 5 kfiJl{w q ) ^ + 2 Attnt + N nntf) cos (qn) 

~ n—1 

+ A (q) (17) 


the Hamiltonian ([!]) of a chain of N = 16 spins with pe¬ 
riodic boundary conditions. Fig. |T| shows the result as a 
function of frequency to for a fixed value of momentum 
q = 7 r along the chain at different temperatures. The 
magnetic fields are h x = 0.1 and h z = 0 . 2 . 

At low temperatures the only visible peak occurs at 
to = J z + h z . It can be attributed to the flipping of a 
single spin in an otherwise all spin-up background. The 
intensity of this spin-wave excitation diminishes as the 
temperature is increased, and new low-energy peaks ap¬ 
pear at w = kh z , where k is an integer. This is the 
finite temperature signature of inter-level WZL transi¬ 
tions. At high temperatures the intensity of the k = 1 
peak becomes the same order of magnitude as the spin- 
wave peak. 



(0 


FIG. 1: (Color online) Dynamical structure factor S xx (q = 
7r, to) obtained by exact diagonalization on a 16-site chain for 
different temperatures indicated by the legends. h z = 0.2 
and hx =0.1. The delta functions have been broadened using 
Gaussians of width A = 0.005. The inset shows the integrated 
intensity of the peaks, S^ x (q), vs q for k = 1 , 0,2 (top to 
bottom) at T = 0.5. 


where w q = 2xsu\(q/2) and A (q) ~ 0(x 2 ) is a correction 
due to extending the sums to negative infinity 

x 2 

A(g) = 32 + N^^) 

w 2 

~ (2A-J-J.-I- + 6Aj-|.^ + -ZV44444) 

x 2 

~ (7A tw + 2Am!!) 

9 

+ / ( 5 A- 14.14 + 21 Vpj-!-|-|-) + 0(x A ). (18) 


V. NUMERICAL RESULTS 


In order to validate Eq.(17) we will calculate the dy¬ 
namical structure factor S xx (q, to) numerically. To calcu¬ 
late S xx {q 1 to) at finite temperatures we have diagonalized 


Some of the main features of the finite temperature dy¬ 
namical structure factor in Fig. [l] can be explained qual¬ 
itatively using what one expects in the limit of vanishing 
h x when the Hamiltonian ([l]) becomes purely classical. In 
that limit the peak at to = h z is associated with flipping a 
spin at a domain-wall. This process does not create new 
domain-walls, and the energy cost is therefore only the 
cost of one extra spin opposing the magnetic held. The 
peak at to = J z — h z can be attributed to hipping a spin 
inside a domain of spins that are pointing opposite to the 
magnetic held. This creates two new domain-walls and 
makes one less spin pointing opposite to the held. These 
two processes can only occur at hnite temperatures as 
they require the presence of domain-walls. Other fea¬ 
tures, such as the existence of peaks at to = 0,2 h z ,J z 
cannot be explained from considerations at h x = 0 . 

When h x > 0 one expects generally that the location 
of the peaks in S xx {q 1 ui ) will acquire a dependence on 
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the momentum q. This is true for the spin-wave peaks, 
but not for the low-energy peaks associated with the 
inter-level WZL transitions. They remain at frequencies 
c j = kh z for all q. However, according to Eq. (17) their 
spectral weights should depend non-trivially on the mo¬ 
mentum q. To estimate numerically the spectral weight 
of each low-energy peak at ui = kh z , k = 0, 1 , 2 , we in¬ 
tegrate the dynamical structure factor numerically over 
a small frequency interval chosen big enough to capture 
the full peak. We choose w £ [(k — l/2)h z , [k + 1/2 )h z \, 
and label the resulting integrated peak intensity S% x (q) 
as in the analytical calculation. A plot of these as func¬ 
tions of q is shown in the inset of Fig. [l] for T = 0.5. 
The k = 1 peak has the biggest spectral weight which 
decreases with increasing q , while the spectral weights of 
the k = 0, 2 peaks increase with q 1 except for the Bragg 
contribution at ui = 0 , q = 0 which is due to a finite 
ground state magnetization in the x-direction and is un¬ 
related to BO. Our analysis is restricted to the lowest 
values of k because at higher values of k there is an over¬ 
lap with the peak at u> = J z — h z ■ To keep these separate 
we require J z — h z > (k + 1/2 )h z . 


VI. COMPARISONS 

To compare the analytical result Eq. with our nu¬ 
merical results we consider each low-frequency peak (k) 
of the dynamical structure factor separately. 


A. k=l 

The expression for S xx (q) simplifies considerably at 

q = 0 , 


Sfc=i(0) — 2^ttl — ^7 ( 71 V.i~i.pr + 2jVj.-j-4.-1~t.) ■ ( 19 ) 

The upper panel of Fig. [2] shows how the expression 
Eq. (19) compares to the exact diagonalization results 


for Si\q = 0). It is rather accurate for x < 1, but un¬ 
derestimates the numerical result for x = 1. According 
to the analytical result the dynamical structure factor for 
k = 1 depends on q through the variable w q . Expanding 
for small w q we find 


SUM = sn 0) - ^ (&N m - 5jV tm - 2jv mtt ) 

+ 0 (x 4 ) ( 20 ) 

The lower panel of Fig. [2] shows the exact diagonalization 
results for Sf x (q) — S xx (0) plotted as a function of w q 
for different values of x and q at the same temperature. 
The data points for x < 1 fall on the same curve, and an 
asymptotic linear dependence on w q , shown as the dotted 
line, is seen for small w q . In the inset we have plotted the 
negative slope of this dotted line, the point at T = 0.5. 



FIG. 2: (Color online) fc = 1 peak for h z = 0.2. The upper 
panel shows exact diagonalization results (symbols) for the q = 
0 intensity vs. T for different values of x = 2 h x /h z . The solid 
curves are gotten by Eq. (191. The lower panel shows how 
the q dependence of the k = 1 peak scales as w q for a fixed 
T = 0.5. The inset shows the negative slope of such scaling 
curves vs. T for x = 0.2. The solid curve is the expression 
(6A4.-J.4. — 5 A mt — 2jV-|~|-4,-|~|-)/iV. 


Repeating the scaling plot for several temperatures we 
get the remaining points. The solid curve shows the an¬ 
alytical result for the slopes, Eq. (20). For T < 0.5 the 
analytical results agree well with the numerical results, 
while at higher T the analytical result overestimates the 
slope. The reason for this is probably that interactions 
between spin-clusters becomes dominant at high temper¬ 
atures and densities and must be treated more accurately. 


B. k=2 


The q = 0 part of the k = 2 peak is small and of 0(x 4 ), 
whereas 5f x ((?) depends on q through the variable w q . 
The leading behavior for small w q is 


s x k iM = ST{ o) + ^i 


W n 




+ 0 (x 4 ). ( 21 ) 


In Fig. [3] we have plotted S% x (q) — S , 2 X ( 0 ) from exact di¬ 
agonalization vs. w q for several values of q and magnetic 
fields at a temperature T = 0.5. The points fall on the 
same curve demonstrating scaling with w q . The dotted 
line shows the asymptotic linear behavior of the curve. 
It was obtained by lumping all the data points into one 
set, and fitting them to a quadratic function. In the in¬ 
set we show the slope of this curve together with slopes 
of scaling functions obtained at other temperatures and 
compare them with Eqs. (16) and (21). As in the k = 1 
case the comparison is good up to T ~ 0.5. 
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FIG. 3: (Color online) Scaling of the spectral weight of the 
k = 2 peak with Wq. The main panel shows 8(Sf :c ((?) — 
Sf®(0)) /N from exact diagonalizations vs. w 2 at T = 0.5 
for several values of magnetic fields indicated by the legends. 
The dotted line is the asymptotic linear behavior for small w 2 . 
The inset shows the asymptotic slopes of the scaling curves 
for h x = 0.02, h z = 0.2 vs. T. For comparison N-ffi/N from 
Eq. (161 is plotted as the solid curve. 


C. k=0 


from intra ground state transitions which cause the dis¬ 
continuous behavior at q = 0 seen in the inset of Fig. [I] for 
k = 0. This additional term is proportional to the ground 
state transverse magnetization squared. Although it is 
the dominant term, we will ignore it here as it is not re¬ 
lated to BO and only occurs at q = 0 and ui = 0. Fig. [4] 
shows a comparison of the numerical results with the ex¬ 
pression Eq. (23) for several values of the temperature 
and ratios of magnetic fields. The agreement is reason¬ 
able for low values and temperatures T <~ 0.5. 



For the zero energy 
the sum over n in Eq. 
sions, Eq. (fl4|), we get 


peak, k = 0, we need to evaluate 
Using the classical expres- 


'y ' (A^-nt "k + Apt-ntt) cos (qn) 

n—1 




cos(g) — K_ 

1 + K' l _ — 2AT cos (q) 

( 22 ) 


Expanding to order 0{x 2 ) we find 
x 2 . 

So x (q) = + 8Aft4,t + ^V’tTFTt) 

+ -g - (Nn + 3iV m ) 

9 

W q 

g- (^ttt “k 2iVt.fj,f + -^t'TTTT) 
cos (q) — I\- 
X l+ K 2 _- 2AT cos {q) 
w 2 

21 (2A r |j_| + + fV|-(-|-|'|). (23) 

This expression is a smooth function of q. Thus there are 
zero energy excitations also at a finite q. In addition to 
this expression, at q = 0, there is a term coming from the 
finite ground state magnetization in the ^-direction, i.e. 


FIG. 4: (Color online) Spectral weight of the k = 0 peak as 
functions of momentum for different values of temperature 
and magnetic field ratios x as indicated by the legends. h z = 
0.2. The solid dots are exact diagonalization results. The 
solid curves are the analytical result Eq. (231. 


VII. COBALT NIOBATE 

The material CoNb206 is an experimental realization 
of a system of weakly coupled one-dimensional Ising-like 
ferromagnets. It has been studied with neutron scatter¬ 
ing and clear evidence of dispersive spin-cluster states 
have been observed in zero applied magnetic field^ 6 . 
The fact that the excitations are dispersive at zero fields 
means that they cannot be interpreted as belonging to 
the WZL. Furthermore, it means that the Hamiltonian 
of CoNboOg must contain extra terms in addition to 
the dominant Ising coupling 16 . A more detailed micro¬ 
scopic model which includes both internal magnetic fields 
and extra spin-spin couplings has been proposed for this 
material 1 '. In particular, the large ratio of the internal 
fields h l x /h l z ~ 10 causes the spin-clusters to be disper¬ 
sive over a large region of momentum space at low ener¬ 
gies, and the WZL-ladder is present only near the zone 
boundarjfl 1 which in the neutron scattering experiments 
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is dominated by the “kinetic bound state” caused by the 
extra spin-spin couplings. 

Still, one should expect to see the finite temperature 
signatures of MBO in CoNb206 if the magnetic fields can 
be arranged so that h x /h z < 1. To investigate this closer 
we have repeated the exact diagonalization of the skew 
field Ising model including also additional termJ^. The 
extra terms we have added to the Hamiltonian ([I]) are 

H xCo m = -Jr E, (S?Sf +1 + S?s? + 1 ) + Jnn E, SfSU 2 

where J x = 0.214 and Jnn = 0.247. The magnetic 
field components, including both internal and external 
magnetic field contributions, have been set to h x = 0.1, 
h z = 0.2, and the temperature is T = 0.5. Fig. [5] shows 
the resulting dynamical structure factor as a function of 
frequency for eight different momenta q £ [0,7r] plotted 
on top of each other. 

At high frequencies, (w 0.6), there are dispersive 
peaks. For low frequencies one can clearly see the two 
non-dispersive peaks at to = kh Zl k = 0,1 associated with 
transitions between WZL states even in the presence of 
the additional couplings. The spectral weights of these, 
and for a smaller peak at k = 2, are shown in the inset as 
a function of momentum. The general trend is the same 
as without the extra couplings, except that the intensity 
of the k = 1 peak shows a larger dependence on q and 
the k = 0,2 peaks have in general a larger intensity. 


VIII. CONCLUSIONS 

We have investigated the finite temperature spectral 
signal of MBO in one-dimensional Ising models in a 
skew magnetic field using both analytical calculations 
and numerical exact diagonalizations. The main neu¬ 
tron scattering signature of MBO are finite temperature 
low-frequency peaks at o o = kh z , where the k = 1 peak 
generally has the biggest spectral weight. While the peak 
frequency does not depend on momentum, its spectral 
weight does. We have shown that the spectral weight of 
the peaks can be calculated analytically when h x /h z < 1 
and h x < T <~ J z /2 using an approach that treats the 
excitations of the chain as a gas of spin-cluster excita¬ 
tions. 

We have also investigated numerically how these re¬ 
sults change when adding other couplings relevant for 
the material CoM^Og, and conclude that the extra cou¬ 
plings present will not destroy the MBO signatures. 

This research was supported by a grant NFR-213606 
from the Research Council of Norway, and the numer¬ 
ical calculations were performed on the Abel computer 
cluster with a Notur grant nn4563k. 


Appendix A: Bessel sums 




CO 


FIG. 5: (Color online) Exact diagonalization results on a 16- 
site chain for parameters relevant for the material CoNb 2 06- 
The main panel shows S xx (q,uj) vs. u> . Each curve is for a 
specific q- value, and q takes values from 0(red) to 7r(violet). 
The magnetic fields are h x = 0.1, h z = 0.2 and the tempera¬ 
ture is T = 0.5. The curves are obtained by broadening the 
delta-functions using Gaussians of width A = 0.005. The in¬ 
set shows the spectral weights S% x (q) for the low-frequency 
peaks with k = 1,0, 2 indicated by the labels. 


The expressions for the matrix elements needed in the 
calculation of the dynamical structure factor involve an 
infinite sum over l. Substituting in the Bessel function 
forms for if), Eq. Q, for x < 1, and omitting the normal¬ 
ization constant which is very close to unity we get 

oo 

A np = ^ E Jl+l-n-k(z p+q )Jl-n^py (ql ~ P)/2 , (Al) 
1=1 
OO 

B np = - £ •/,_!_„_*(**.,)(A2) 
^ 1=2 


In order to evaluate these sums we add and subtract 
terms with negative values of l. For A np we get 

1 OO 

A n p = ^ ^ ^ Jl+l—n—k(Zp-\-q) Jl—n^Zp)^ ^ ^^ ^A n p 

l =—oo 

= ^ J 1 - k (w q )e i ^ n+ < 1 - k ^/ 2 + A A np (A3) 

where w q = 2xsin(g/2) and the sum has been evaluated 
using Graf’s summation formula 1 ^, and A A np is correct¬ 
ing for the added terms with l < 0: 


= E Jl+l-n-k(Zp+l)Jl-n(Zp)e iiql - p)/2 

/=—OO 


<W(-l) fe 


Z P+q Z P -ip /2 

2 k + 2 k\ 


(A4) 














where we have replaced the sum by its lowest order term 
in x which is 0(x 2n+k ~ 1 ). In a similar way 


B np = -- J_ 1 _ k (w q )e« ,in+ « 1 -Q- kp V 2 + A B np , 


(AS) 


v i+fc 
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